Method to estimate ISCO and ring-down 
frequencies in binary systems and consequences for 
gravitational wave data analysis 

Chad Hanna 1 , Miguel Megevand 1 , Evan Ochsner 2 , and Carlos 
Palenzuela 1 ' 3 

1 Department of Physics and Astronomy, Louisiana State University Baton Rouge, 
LA 70802, USA 

2 Department of Physics, University of Maryland, College Park, MD 20742 

3 Max-Planck-Institut fur Gravitationsphysik, Albert-Einstein-Institut, Golm, 
Germany 

E-mail: channa@phys.lsu.edu, megevand@phys.lsu.edu, evano@umd.edu, 
carlos@lsu.edu 



Abstract. 

Recent advances in the description of compact binary systems have produced 
gravitational waveforms that include inspiral, merger and ring-down phases. 
Comparing results from numerical simulations with those of post-Newtonian (PN), 
and related, expansions has provided motivation for employing PN waveforms in 
near merger epochs when searching for gravitational waves and has encouraged the 
development of analytic fits to full numerical waveforms. The models and simulations 
do not yet cover the full binary coalescence parameter space. For these yet un-simulated 
regions, data analysts can still conduct separate inspiral, merger and ring-down 
searches. Improved knowledge about the end of the inspiral phase, the beginning of the 
merger, and the ring-down frequencies could increase the efficiency of both coherent 
inspiral-merger-ring-down (IMR) searches and searches over each phase separately. 
Insight can be gained for all three cases through a recently presented theoretical 
calculation, which, corroborated by the numerical results, provides an implicit formula 
for the final spin of the merged black holes, accurate to within 10% over a large 
parameter space. Knowledge of the final spin allows one to predict the end of the 
inspiral phase and the quasinormal mode ring-down frequencies, and in turn provides 
information about the bandwidth and duration of the merger. In this work we will 
discuss a few of the implications of this calculation for data analysis. 
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1. Introduction 

Gravitational radiation emitted during the inspiral, merger and ring-down of two 
compact objects with a total mass below ~ 1OOM is a likely detection source for 
interferometric gravitational wave detectors pQ. Maximal signal-to- noise ratio (SNR) 
is achieved in searches for compact binary coalescence with matched filtering banks of 
template waveforms [2j |3j HJ 13 E] allowing faint signals to be detected reliably. The 
SNR depends on having a template waveform that matches the true signal correctly 
[2] . Because templates are weighted by the noise power spectral density of the detector 
with non-trivial frequency dependence it is important when in the binary evolution the 
template matches well. Numerical simulations are improving our knowledge of the full 
waveforms emitted during the inspiral, merger and ring-down phases of binary black 
hole coalescence El |9], yet they still do not cover an adequate parameter space for 
constructing filter banks [101 |9]. When the full parameter space is simulated it will 
be useful to have fully parameterized models of these waveforms and even now many 
have considered hybrid waveforms that patch together PN and numerical simulations 
[TT| [T2l El [El E] as data analysis templates. Others have introduced a phenomenological 
fourth order correction to make effective one-body (EOB) waveforms match numerical 
relativity results near the merger [TT| [7] known as pseudo 4PN (p4PN). 

Before the full parameter space is simulated it is still possible to handle the three 
binary coalescence epochs - inspiral, merger and ring-down - separately. We will show 
that a recent calculation [T3] which has power to predict the final spin of a binary merger, 
can provide both insight into when the inspiral phase ends and what the quasi-normal 
ring-down frequency of the merger product should be for a large, yet un-simulated, 
parameter space. We will discuss the following consequences of this calculation: 

• The calculation of the final spin given in [T3] contains an implicit reference to the 
innermost stable circular orbit (ISCO) of a test particle orbiting the Kerr black 
hole that results from the merger. We will also discuss the impact of using the 
new ISCO frequency as a cut off for current gravitational wave searches that are 
targeting coalescing compact objects including the impact on SNR and the ability 
to probe neutron star properties. 

• Knowing the final spin gives an estimate for the quasi-normal mode ring-down 
frequencies. These frequencies could be used in a separate matched filter search 
[26]. 

• We can also use the final spin to calculate the light ring (smallest, unstable circular 
orbit for a photon orbiting a Kerr black hole) and quasinormal mode ring-down 
frequencies of the final black hole; we will see that the 1=2, m=2, n=0 QNM is quite 
close to the light ring frequency. We will define a merger epoch as the evolution 
from ISCO to light ring and dicuss how to obtain its time frequency content for yet 
unsimulated systems. 

We define the merger epoch as everything that occurs between the ISCO andi ring- 
down frequencies. We show that the duration of the radiation seems to be well 
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described by the infall duration of a test particle orbitting the merged black hole 
just below the ISCO. Knowing the duration and bandwith gives a time-frequency 
volume, which is useful in data analysis for unmodeled signals. 

2. Formalism 

The inspiral phase of two compact objects is generally considered to end when the binary 
has evolved to an innermost stable circular orbit (ISCO), should it exist. The ISCO can 
be defined by waveform models and numerical relativity p2J [18j [7] . Unfortunately there 
is a difference in the numbers obtained from different methods. Some LIGO searches 
for non-spinning systems have taken the conservative approach to use the ISCO defined 
for a test particle orbiting a Schwarzschild black hole. 

Intuitively the space-time of an inspiralling comparable mass binary system will 
not be well described by a Schwarzschild space-time, and thus the Schwarzschild ISCO 
will not capture the dynamics correctly. If the merger product is a black hole, it will be 
a Kerr-like black hole with a final spin that is a function of the masses and spins of the 
components. Recently [13] has proposed a simple set of assumptions that predict the 
final spin of the black hole merger product to well within 10% of numerical simulations 
based on first principle arguments. Implicitly in the derivation of [13] is the ISCO radius 
of a test particle orbiting a Kerr black hole having a final spin equal to the spin of the 
merger product. The fact that the ISCO radius used gives the correct answer for the 
final space-time suggests that it plays an important role in the pre-merger dynamics. 
We hypothesize that the test particle ISCO of the merged Kerr black hole may describe 
a way to define the end of the inspiral phase for binaries that produce black holes. 

The ISCO solution for a test particle orbiting a Kerr black hole is |2Uj . 

nsoo = M{3 + Z2=F[(3-Z,)(3 + Z 1 + 2Z 2 )] 1 / 2 } , (1) 

where M is the total mass of the black hole and a/ is the angular momentum. The 
final spin is required in order to compute this ISCO. In [13] . assuming the amount of 
mass and angular momentum radiated beyond the ISCO is small, the following implicit 
formula for the final angular momentum of a black hole a/ with component spins aligned 
with the orbit is calculated, 

a/ = Lorb ( Z__ r isco o/\ , 1 2 Xi + X2 (2) 
M M 2 v' M M ' MJ (1 + q) 2 ' U 

where x% = di/n^u Q = m,\jmi G [0,1] and M = mi + m 2 is the total mass. The 
implicitly found a/ agrees well with numerical simulations [13] and the analysis can 
be modified to include arbitrary spin angles. L or b is the orbital angular momentum 
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contribution calculated from the orbital angular momentum of a particle at the ISCO 
of a Kerr black hole with spin parameter a,f, which has the following expression [20] . 

iW_ r a,\ g ±(r 2 T 2a / M 1 /V/ 2 + af 



M 2 V y 'M'M7 (l + g) 2 M 1 /2 r 3 / 4 (r 3 / 2 -3Mr 1 /2 ± 2a / M 1 /2)i/2 ' ^ > 

where the upper/lower signs correspond to prograde/retrograde orbits. In order to 
agree with numerical simulations this function has to be evaluated at r = rjsco given 
by equation Q. Alternative radii in the above prescription - like the photon- ring 
radius or Schwarzschild radius - give a prediction that is off the numerically obtained 
value. This fact, combined with the same analysis extended to elliptic binaries, again 
depending sensitively on the ISCO location [23], suggests the ISCO radius determined 
by this approach is relevant [24] . 

In order to get the gravitational wave frequency at ISCO we use the coordinate 
angular velocity of a circular orbit 



fi = ± r3/2± G/ MV2 ' ( 4 ) 

with <if given by the implicit equation ([2]). The gravitational wave frequency at a given 
radius is then / = Q/n, and so we define /isco [bkl] as the frequency obtained by solving 
the system of equations (|l]|4]) at the Kerr ISCO radius. 

The solution space for /isco [bkl] can be written as a function of the final unknown 
spin cif. For convenience, we prefer to extend it as a surface parametrized by (a/,g,x), 
as it is shown in figure [T} This way the mass ratio dependence of the lines corresponding 
to different individual spins can be seen explicitly. For the equal mass case without spin 
the approximate expression for the Kerr test particle ISCO frequency is, 

/isco [bkl] (M, q) « (0.8g 3 - 2.6g 2 + 2.8g + l) x n ^ )3/2 , (5) 

which can be compared to the Schwarzschild test particle ISCO /isco[sch] = 

TT-^M)- 3 / 2 . 

Once the final spin is known it is possible to compute the last stable photon orbit 
("light ring") and the quasi-normal mode ring-down frequencies [16J. The light ring for 
a Kerr black hole is, 1201 



rught = 2M{1 + cos [2/3 cos" 1 (Ta f /M)}} (6) 



r 



±a f M 1 / 2 



light =■= «/ 

An approximate fit to the quasinormal mode ring-down frequencies as a function of the 
final spin is given in [25l [26| 127] . For the 1=2, m=2, n=0 mode we have, 

^qnm « °^ (1 - 0.63(1 - a f ) - 3 ) , (8) 

Others ring-down modes could could be computed as necessary. 

Figure [2] compares the BKL ISCO for non-spinning binaries with the ISCOs 
computed from the minimum energy conditions at 3PN order [TH] as well as 2PN and 
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Figure 1. The surface of solutions of the frequencies at the ISCO as a function of 
the mass ratio q and the final spin a t for components with spins that are aligned with 
the orbital angular momentum. Also shown, are curves corresponding to the solution 
of the equal spin case xi — X2 = X- 

p4PN EOB ISCOs [T51ITS]. BKL and EOB agree exactly with the Schwarzschild ISCO at 
the extreme mass ratio limit, but the PN calculations are known to be inconsistent [IH] 
in that regime. Predicting the final spin gives the / = m = 2 ring-down mode frequency 
by ([8]) and the light ring frequency. Both are shown in figure [2] The ISCO and light ring 
set a natural "merger" epoch which can be searched for by burst search techniques (28] in 
between the inspiral and ring-down matched filter searches. The ring-down frequencies 
shown in [16] agree well with these predictions. Although figure [2] is plotted for the 
non-spinning case the formalism described above can be easily generalized. 

It has been shown j29j [301 E2] to be useful to understand un- modeled time-frequency 
evolutions in gravitational wave analysis in terms of time-frequency volumes, Af At. 
This concept can provide insight into the merger epochs of various systems. We will 
consider three timescales for the merger 1) the Newtonian free fall time scale tg, 2) the 
quadrupole radiation time scale t qua d and 3) the time scale of a test particle following 
its geodesic just inside the ISCO of the Kerr space-time (corresponding to the merger 
product). For comparison we also plot the time- frequency volumes found from the 
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Figure 2. ISCO and QNM frequency estimates for non-spinning binaries as a function 
of mass ratio using different methods. Knowing the final spin of the black hole gives the 
expected ring-down frequency which agrees with the light ring frequency. The ISCO 
frequencies are very different depending on the method. The PN minimum energy 
condition gives an inconsistent result in the extreme mass ratio limit whereas the the 
other methods (EOB, BKL) agree with the Schwarzschild ISCO at small values of q. 
The ISCO and QNM frequencies define a natural merger epoch which can be analyzed, 
albeit sub-optimally, even without knowing the numerical waveform. 



figures in [IB]. The Newtonian free fall time-frequency volume is, 

2 -3/2 

A/ ISCO At ff = . (9) 

y/ivniisco 

The quadrupole time-frequency volume is calculated by assuming that an ISCO doesn't 
exist and that the system continues to be driven by radiation. To calculate this we use 
the quadrupole approximation to the inspiral waveform given by [21], which gives the 
frequency evolution as, 

"5(1 + qf M 13/8 



256 q t- t 
This leads to the time frequency volume, 



(10) 



A/bklisco AVad = 5 ^ fi +9) [(Mft Isco r 8/3 - (Mfi light )- 8 / 3 ] 

x [(Mfi %ht ) - (Mflisco)] , (11) 

where the ISCO frequencies and light ring frequencies are obtained from the BKL 
approximation as in figure [2] 
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The merger evolution will not be purely dominated by radiation by virtue that 
there are no circular orbits below the ISCO. The numbers obtained by this estimate 
should greatly over predict the time scale for extreme mass ratios where particles would 
take an asymptotically infinite time to fall in. The free fall time scale will under predict 
the time since the system may still complete some (unstable) orbits before merging. 
The estimate that agrees best with numerical relativity simulations is the infall time 
of a test particle falling from a circular orbit just below the ISCO to the light ring 
of the the merged Kerr black hole. Figure [3] shows the time-frequency volume for 
the merger epochs as a function of mass ratio for non-spinning binaries. The figure 
shows that the numbers taken from the simulations presented in [IB] agree well with 
the predictions from the test particle in the Kerr space-time. By these estimates the 
mergers from yet un-simulated evolutions (low mass ratio) are likely to have the largest 
time frequency content and searches may benefit from the numerical simulations of these 
[30] . This information should be useful in conducting IMR searches and also in guiding 
the construction of analytic full waveform models. It is of course possible to repeat 
this analysis for arbitrary spin configurations using the above prescription and others 
[T51 12T| [221 E3] • We will leave further discussion of these matters to future work and 
will spend the rest of this work considering the implications of using the ISCO found 
from the above set of equations as a termination frequency for PN waveforms in initial 
LIGO data analysis. 




0.01 0.1 1 

Mass ratio, q 



Figure 3. An estimate of the time-frequency volume for the merger epoch of non- 
spinning binary black holes. Numerical relativity results agree well with the time- 
frequency volume of a test particle falling into the merged Kerr black hole from a 
circular orbit just inside the ISCO. 
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3. Possible impact on data analysis for the inspiral phase of compact 
binary coalescence searches 

As mentioned previously, there is ambiguity on how to define when the inspiral phase 
ends. However, what is really necessary is the characteristic frequency at which a given 
template waveform ceases to resemble the numerical simulations. Some searches for low 
mass binary systems of non-spinning component masses use the ISCO of a test particle 
orbiting a Schwarzschild black hole [Hj. Other searches use the ISCOs calculated from 
explicit PN energy considerations, and some abandon the use of an ISCO altogether 
and use the Schwarzschild light ring as a termination frequency [15J. 

Recall that for low mass ratio (q ~ 0) systems of non-spinning objects, the 
Schwarzschild test particle limit is a good approximation for the expected ISCO 
frequency since the merger product will be a Schwarzschild- like black hole. However, for 
systems with comparable masses and/or spin the true ISCO frequency may be different 
since the non-trivial contribution from the orbital angular momentum will have a strong 
impact on the final black hole's spin, and the space-time in the near merger epoch. Many 
have addressed how well various PN approximants stay faithful to numerical relativity 
solutions[IU HU H2, 03]. Most approximants stay faithful through the Schwarzschild 
ISCO frequency [13]. Some approximants fare extremely well beyond this point [43]. 

The fact that some approximants do remain faithful far beyond the Schwarzschild 
ISCO leads us to examine using the ISCO frequency described in the previous section, 
equation ([5]), as the termination frequency for inspiral data analysis. As shown, this 
frequency is consistent with some of the PN and EOB models for predicting the ISCO in 
the equal mass regime (fig. [2| and is consistent with exact solutions in the test particle 
limit. It also has the advantage of being waveform model, or fit, independent (based 
on first principles) and is easy to calculate. We will conservatively model the impact 
of possible phase errors incurred by using an approximant that doesn't quite match 
numerical relativity results in our assessment. 

We use the stationary phase approximation waveforms, which are often employed 
in searches, and can be defined as, 

h(f) OC /-7/6 e i*(/;M,„) 5 (12) 

where M is the total mass and \i is the reduced mass [H]. In order to ascertain the 
phase fidelity we turn to comparisons given in [11] which show for equal mass binary 
simulations the phase difference between numerical results and 3.5 PN waveforms as 
a function of frequency. The figures in [TT] indicate a good fidelity in phase through 
the Schwarzschild ISCO frequency. However substantial phase errors may accumulate 
between the Schwarzschild and BKL ISCO frequencies. In order to accurately predict 
the impact of integrating PN waveforms to the BKL ISCO we decided to model 
the possibility that secular phase errors could be ±7r or more radians between the 
Schwarzschild and BKL ISCOs. We propose calculating the SNR ratio of an ensemble 
of waveforms where phase errors are allowed to accumulate linearly between the 
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Schwarzschild ISCO / 2 (M,g) and the BKL ISCO f x {M,q) given by, 

snr /2 " \ U/ low f^SM) ) W /low Wsjf) 

where R n is a normally distributed variable with mean and variance 1 . fi ow is 40Hz [5] . 
[|] We note that the above calculation is a conservative estimate of the SNR gain for two 
reasons: 1) A template with different parameters may match the numerical waveforms 
better, producing better phase agreement. Since search results are maximized over SNR 
this is an important possibility. 2) The phase errors do not accumulate linearly between 
these frequencies, most of the phase error occurs near the BKL ISCO. Using the {M, q] 
dependent ISCO frequencies as fi(M,q) and the Schwarzschild ISCO frequencies as 
/ 2 (M) we plot the expected SNR ratio for initial LIGO [35] in figure |1J which shows that 
for some total mass and mass ratio combinations there is an appreciable gain in SNR by 
integrating to the new ISCO frequencies despite modeled phase errors. This calculation 
is not intended to suggest that SPA templates cut at the BKL ISCO produce the optimal 
gain in SNR, since other approximants or frequency cut offs could do better. Instead, 
it aims to simply illustrate the effect of integrating to the BKL ISCO as a function of 
M and q for present searches. The mass range is limited to 80 M & so that at least two 
cycles exist in the waveform between 40Hz and the Schwarzschild ISCO frequency. 




4. Probing the tidal disruption of neutron stars 

If the BKL ISCO truly marks the transition from an inspiral, radiation-dominated 
evolution to a dynamical one then it does have some consquence for the ability to 
probe the tidal disruption of neutron stars through gravitational waves. Estimation of 
gravitational-wave frequencies for tidal disruption of NS-BH binaries may be useful for 
determining various properties of the neutron star, such as the radius (36], |37] . Knowing 
the radius of the neutron star would in turn provide information about its equation 
of state [381 ES]. However, as pointed out in |36j, it would be difficult to extract 
information about the disruption unless it occurs before the binary reaches the ISCO. 
If the disruption occurs afterwards, its signature in the produced gravitational waves 
might be too weak for extracting accurate information. Thus, the "useful" cases are 
those for which the tidal-disruption frequency is less than the frequency at the plunge 

{ftd < /plunge = /iSCo)- 

A more accurate calculation of f p i unge would imply a better estimate of the range 
of NS-BH binaries for which the disruption is more plausible to be measured. In this 
section we show how that range changes when the formalism of section [2] is applied, as 
compared to that in [36] . 



| We take S n (f) to be approximated analytically as 



S n (f) « 9 x 10 



-46 



~ U!)/ ' ' 0.1q(JJ) +0.32(J-\ +0.52 



150 7 V 150 / V 150 



(14) 
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Figure 4. The SNR ratio defined by equation ( [13) for the non-spinning case x = 0. 
The solid lines include phase errors of the form ( |13| whereas the dashed lines assume 
perfect phase coherence for comparison. 



The frequencies obtained in [36J correspond to the extreme mass ratio case, q = 0, 
and non-spinning neutron star, %ns = 0. In figure [5] we reproduce those results 
(compare to figure 2 of [US]), together with the values obtained using the formalism 
presented in this work. We also assume in this case that the neutron star is non- 
spinning, but with q = m NS /m B H 7^ 0. The curves of R vs / t( j are calculated using 
the formula given in equation. (3) of [36], for the same parameters as in that work: 
mNs = 1.4M ;mBH — 2.5, 10, 20, 40, 8OM . The circles over each curve specify where 
the plunge occurs for different values of Xbh- Filled circles correspond to the extreme 
mass ratio approximation. The meaningful part of each curve will be that to the left of 
the white circles (f td < / p i un ge)- 

The plunge frequencies calculated in |36j for the case with tubh = 2.5M© (black 
circles in the right-hand-side curve of figure |5| correspond to retrograde orbits. However, 
for that mass ratio (q = 0.56), the test particle orbiting the merger product ISCO does 
so in a prograde fashion due to the orbital angular momentum dominating the final 
spin. Thus, the final spin is aligned with the orbital angular momentum even though 
its initial spin is anti-aligned. 

We now compare with the NR results obtained in [40J. We calculate /bkl for the 
six cases considered in in that work, models A-F. The BKL frequencies are shown in 
table [TJ together with the corresponding parameters q and M = Mbh + -^ns (xbh and 



11 




400 600 800 1,000 1,200 1,400 1,600 1,800 2,000 
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Figure 5. Radius of the neutron star vs. the disruption frequency. The circles 
indicate the points at which the plunge would occur; the filled circles correspond to 
/plunge as calculated in [35], included here for comparing. The value of xbh is indicated 
next to the corresponding circle where possible, otherwise a dashed line was used to 
connect the value to the circle. Each curve is labeled with the value of the BH mass: 
tobh = 2.5 - 80M o . 

Xns are set to zero as in [ID])- Although the models studied in [10] expand only a small 
range of parameters; we see that, at least in that range, the frequencies obtained with 
the BKL approach are consistent with the NR results (See figure 7 of [10]) - that is they 
predict in every case a higher value than the measured / cut caused by the tidal effects. 

Table 1. BKL frequencies for the same parameters as models A-F of [3D]. These 
values are consistent with the NR results - that is the BKL ISCO frequency is higher 
than the observable / cut from the simulations |40j 



Model q M(M ) / B K L (kHz) / cut from gD] 



A 


0.327 


5.277 


1.38 


1.16 


B 


0.327 


5.244 


1.39 


1.41 


C 


0.328 


5.311 


1.37 


0.92 


D 


0.392 


4.623 


1.65 


1.14 


E 


0.392 


4.594 


1.66 


1.40 


F 


0.281 


5.929 


1.18 


1.09 
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5. Conclusion 

The estimate for the spin of the black hole that results from BH-BH, BH-NS, and 
some NS-NS mergers found in [13] leads to an estimate for the ISCO, light ring, and 
quasinormal mode ring-down frequencies of a general compact binary system that is 
waveform model independent. Both frequencies have impact on searches for the inspiral, 
merger and ring-down epochs of the gravitational waves emitted by these systems. We 
have shown that the most interesting merger-epoch time-frequency volumes are for 
extreme mass ratios. More work may have to be done to model merger epochs for 
these systems. The formalism presented here to describe inspiral, merger and ring- 
down epochs can be extended to arbitrary spin configurations and may help guide the 
construction of better analytic models. 

As for inspiral epoch searches alone, we have shown that for symmetric mass, 
non-spinning systems with greater than 45 M Q total mass there can be a 30% gain in 
SNR over using the Schwarzschild ISCO radius as is done currently [14] even including 
conservative phase errors. It is worth noticing that more general fittings can be 
computed in order to cover the full range of parameters (q, X11X2) within the BKL 
approximation. Additionally we discussed cases where neutron star tidal disruption 
may be observed more easily through gravitational waves. We would like to stress that 
this work is open ended. Here we have only considered spins that are aligned and anti- 
aligned with the orbital angular momentum. It is certainly possible to extend this work 
to make predictions about more general cases, all of which should be compared with 
numerical relativity results as they are produced. 
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